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Abstract 

We analyze the solution of the coined quantum walk on a line. First, 
we derive the full solution, for arbitrary unitary transformations, by using 
a new approach based on the four "walk fields" which we show determine 
the dynamics. The particular way of deriving the solution allows a rigor- 
ous derivation of a long wavelength approximation. This long wavelength 
approximation is useful as it provides an approximate analytical expres- 
sion that captures the basics of the quantum walk and allows us to gain 
insight into the physics of the process. 



1 Introduction 

Consider a two-state particle, e.g. a two-level atom, moving along a one- 
dimensional lattice. The displacement is conditional; that is, the particle jumps 
one site at each time to the right or to the left depending on its internal state. 
The movement of such a particle is rather trivial: If the particle is in its "up" 
or "down" internal state it moves to the right or to the left, respectively, and it 
undergoes a superposition of right and left displacements if it is in a superpo- 
sition of its "up" and "down" states. But the movement becomes nontrivial if 
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after each displacement a suitable unitary transformation, e.g. the Hadamard 
transformation, leads to a superposition of the internal states of the particle. 
Now the probability distribution that the particle would be found at a given po- 
sition is quite rich: It is null at odd (even) positions at odd (even) times, shows 
two strong peaks near the right and left edges, develops strong oscillations near 
these peaks and has a smooth plateau in the central region. This process is 
known as the (coined) quantum walk. 

In a sense, the quantum walk, which we refer to as QW, is the quantum 
generalization of the classical random walk, and as such was first introduced 
ten years ago by Aharonov et al. Pp. It was later rediscovered by Meyer j2| in 
the context of quantum cellular automata, and more recently by Watrous |3] 
in the context of quantum algorithms. The QW introduced in these papers is 
called the discrete QW (be it coined or not [3]), and must be distinguished 
from the so-called continuous QW, introduced by Farhi and Gutmann in 1998 

which realizes to some extent a quantum version of continuous time classical 
Markov chains. The connection between these two types of QWs is still unclear 

The QW, both in its discrete and continuous versions, is a process that has 
recently attracted the attention of many workers in the quantum computing 
community; see, for example, the recent review by Julia Kempe [5]. Of interest 
is the possibility that QWs will give rise to new quantum algorithms that show 
clear speedups over classical algorithms. The first results in this direction have 
been obtained by Shenvi et al. [6j, who have shown that a QW can perform 
the same tasks as Grover's search algorithm 0, and by Childs et al. |H], who 
have introduced an algorithm for crossing a special type of graph exponentially 
faster than can be done classically. 

In the present communication we give an alternate description of the discrete 
QW that we feel gives physical insight into the nature of that process. This 
paper is an extention of our previous work j9j , in which we showed that the QW 
on the line is not a characteristically quantum process in the specific sense that it 
can be understood simply as a wave phenomenon, thus opening the possibility of 
entirely classical implementations jJOl- Indeed, wave implementations fill I12| 
of the QW have already been discussed. In other words, the characteristic 
quantum feature of entanglement does not seem to play any central role in the 
QW, at least in the unidimensional case. The situation could be different for 
higher dimensional QWs; see the discussion in Knight et al. ^Hl for some insight 
into the role of entanglement in higher dimensional QWs. 

The rest of the article is organized as follows. In Section 2 we give the 
evolution equations for the standard discrete coined QW, and reformulate them 
in a way that directly exhibits an important property: The two coin states evolve 
independently of each other, the coupling between them entering only through 
the first iteration of the walk. In Section 3 we derive the formal solution of the 
Q W in a manner different to those already carried out ^| . This alternative 
derivation allows a rigorous derivation of the long wavelength approximation 
that corresponds to the continuous limit of the QW we pointed out earlier 
In Section 4 we give the solution of this long wavelength limit and discuss some 
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of its properties. Finally in Section 5 we give our conclusions. 



2 Equations of evolution 

Let Hp he the Hilbert space of the particle positions and {|m) , m S Z} a basis 
of Hp; and let He be the Hilbert space describing the two states of the qubit 
associated with the internal states of the particle (usually referred to as the 
coin), and , IL)} a basis of He- The state of the total system belongs to 
the space H = He ® Hp. The dynamics of the system is governed by two 
operations, the conditional diplacement 

D\m,R) ^ \m+l,R), (1) 
D\m,L) = \m-l,L), (2) 

and the transformation acting on the internal states of the coin 



H\m,R) = ^\m,R) + ^n^\m,L), (3) 
H\m,L) = ^/T^\m,R)-^\m,L), (4) 

< p < 1, which is the most general unitary transformation that one needs to 
consider ^Sj (notice that H is the Hadamard transformation when p = 1/2). 
The QW is implemented by the repeated action of the operator HD, i.e. if 
denotes the state of the system after n iterations then 

|V)„ = (^£))"|^)o, (5) 

which can be written as 

+n 

IV')„ = 51 [Rni,n\ni,R) + Ljn^n\m,L)], (6) 
m— — n 

with the equations of evolution of the probability amplitudes given by 



Rm,n — \fpRm-l,n-l + V 1 ^ P-^m+l,n-lj (7) 
im,n = — pRm-l,n-l — \/p^m+l,n-l- (8) 

Finally, the probability of finding the particle at position m at iteration (time) 
n, is given by 

Pm{n) = Pr'l{n)+P.^in) (9) 
P«(n) = \Rm,n\\ P,i{n)^\L, 



^ni,n I 



Eqs.JT)) and © are the discrete QW equations as can be found e.g. in the 
work of Carteret [Tl! for the specific case of the fair coin p = 1/2. Different 
equations of evolution, but still leading to the same final probability, can be 
found by defining the walk as the repeated action of the operator DH, instead 
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of HD, as is done by Nayak and Vishwanath ^3]. In the Appendix we identify 
the transformations that connect these different conventions for the probability 
amplitudes. 

There is a crucial feature in Eqs.Q and © that is at the heart of our 
treatment of the QW. Substituting Eq. © into Eq. (Tj) and rearranging the terms 
one easily gets 



7n-\-2.n — 2 • (10) 

Noticing that the last term in the right-hand side is nothing but -Rm+ i.n-i, and 
making the change n n + 1, we can write 

Rm,n+1 = Rm,n-1 + y/pRm-l,n — \/7>Rm+l,n- (H) 

A similar equation for L^^.n can be derived. That is, Eqs.Q and © are equiv- 
alent to 

Clm,n+1 — am,n-l = ^/p[a7n-l,n ~ am+l,n] , a = R, L (12) 

which is a remarkable equation, as it shows that the evolutions of the two coin 
states are independent. That is, Rm.n does not depend on Lm,n at previous 
times, and vice versa. The only coupling between Rm,n and Ljn,n appears at 
the first interation. Notice that in Eq. (|12l) one needs to specify two initial 
conditions (both 0^,0 and a,„_i) and one needs Eqs.Q and (jS)) for calculating 
am,i from a^.o- But apart from this initial coupling, the evolution is completely 
independent for the two coin states. In other words: The QW is composed of 
two independent QWs, one for each coin state, coupled only at initialization. 

Eq. 1)12(1 is interesting for another reason already discussed in 9 in the special 
case of the Hadamard transformation (p — 1/2). Simple inspection leads to 
the idea that Ea. (|12|) can be understood as the discretization of a differential 
equation involving all the odd derivatives with respect to space and time of a 
continuous field a (x, t) 

^ (2fc+l)!at2fc+i"^'''^^ " (2fc + l)!9x2fe+i"^'''*^- 

This is too nai've an approximation, as the above equation does not respect 
the symetries of the solutions of Ea. (|12|l : Ea. (|13|l describes the propagation of 
waves in the positive position x direction, whilst the probability distribution of 
the quantum walk clearly shows propagation in both the positive and negative 
directions. But it is clear that Ea. (|12|l suggests the existence of some kind of 
continuous description of the discrete coined QW. Earlier ^ we obtained such 
a continuous description by introducing, without justification, new amplitudes 
„ through the definition 
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and then we derived continuous differential equations for the continuous fields 
{x,t). These equations are similar to Ea. (|13|l but with a velocity that differs 
in its sign for the fields A+ {x, t) and {x,t). From these equations we obtained 
an approximate evolution equation by neglecting higher order derivatives 
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A^{^,t), (15) 



with ^ = x/Ax and r = t/At. Eq.JTSl can be solved explicitly in terms of Airy 
functions, and shows that to this degree of approximation the quantum walk 
fields A^ „ satisfy the same equation as the evolution of an optical pulse in a 
medium with third order dispersion. 

The weak points of this approach [5| are that the fields A^ are introduced 
without justification, that only the special case of the Hadamard transformation 
is considered, and that the continuum limit is not clearly derived. The goal of 
this paper is to present a rigorous derivation of Eq. H13|l for arbitrary values of 
p. For that purpose, we first turn to a new derivation of the full solution of the 
quantum walk on the line. 



3 The exact solution: an approach based on four 
fields 

As at a given time identified by the index n the amplitudes am,n depend on the 
discrete variable m. So it is usual jl3[ll4| to define a continuous function of k, 
periodic over a range of 27r, such that the am,n are the Fourier components of 
that function fn{K), 

oo 

U{k)= J2 «m,ne^™"- (16) 
m— ~oo 

It is then possible to formulate the problem in terms of the continuous function 
/„(k) of K rather than the discrete function am.n of m. A similar strategy can be 
implemented with respect to the dependence on the time step n. We essentially 
follow this approach below, but our perspective is slightly different. Because we 
want to be able to understand Eq. (|12() as the discretization of certain linear 
partial differential equations, we introduce Fourier expansions such as in a 
more physical way. 

3.1 Fourier amplitudes 

From the amplitudes am,n defined at discrete points in space, labeled by m, and 
time, labeled by n, we can formally introduce a field a{x, t) at arbitrary points 
in space and time defined by 

oo 

a(x,t)= am,n5 {x — mX) 5 {t — nT) , (17) 

m.n— — oo 
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which produces appropriately weighted delta functions at the positions in space 
and time where Om.n is defined. Here X is the distance between the sites in the 
one-dimensional lattice, and T is the time between steps in the walk. We can 
formally introduce the Fourier transform of this function 

1 /■+°° - - - 
a{x,t)^-—-r dkdu)a{k,uj)e'^'=e-"^\ (18) 

where from Eg. ljl?!) we find that 

/ + 00 _ +00 
a(a;, i)e-*'=^e*^* = ^ a™,„e-^*='"^e*""^. (19) 
-OO ™ ^_ , , 



m.n— — oo 



Next, note that we can write each k and Cj uniquely as 



~k = k + ^, (i = c^ + ^, (20) 
with p and q integers, and k and lo lying within the ranges 

Then using the identity 

+ 00 +00 

J2 S{y-27Tp) = — ^ e^Py, (22) 

P— — 00 P— — 00 

which follows immediately from identifying the Fourier sum representing the 
periodic function of y on the left-hand side, we easily find from Ea. p9|) that 
a{k,Lo) — a{k,uj). Then, making use of Eas. H17|l - H19ll and (|22|l . we can write 

OO 

a{x,t) ^d{x,t) S{x-mX)S{t-nT) (23) 

where 



m,n— — OO 



1 p+tt/X f+T^/T 

a{x,t) = ^XT / dk duja{k,uj)e'''''e~'"^* (24) 

OO 

a{k,u;) = J2 (^rn.ne-"e'^''^, (25) 

m,n— — OO 

cf. Eq.dUl). 

3.2 Dispersion relation 

Returning now to the discrete QW Eas. (|12() . we multiply each side by exp {iujnT) 
exp {—ikmX) and sum over m and n. This yields the equation 

sinter = ^ sin fcX, (26) 
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which identifies the dispersion relations lo — uj{k). Because of the range of 
k and uj appearing in Ea. (|24() . we need only consider such solutions over the 
ranges (I21|l . There are then two such solutions. One solution, which we label 
a;o(fc), is that for which ^ as A: ^ 0. The second, which we label wi(fc), is 
given by 

^i(fc) = -^o(fc)±^, (27) 
for k respectively > and < 0. These relations are represented in Fig.l. 

3.3 Separation into two fields 

Solutions of Eas. ((T^ must respect the dispersion relations, and thus for Eq.(|23J) 

to represent an actual solution we must have 

a {k, Lo) = [a+ (fc) 5{lo~lo^ [k)) + a~ (k) 5{lo-uj^ (fc))] , (28) 

for some amplitudes a+(fc) and a^(k). Thus from Ea. (|24|l we find 

r+Tr/X p+tt/X 

a{x,t)^ dfca+(fc)e*'=="e"*'^°('=)*+ / dfc a" (fc) e^'^^e-^'^^^'^)*. (29) 

J-tt/X J-tt/X 

Using this expression for a{x, t) in Ea. H23|l for a{x, t), employing Ea. H27|l to write 
wi(fc) in terms of uJo{k), and recalling that exp (±i7rn) = (—1)" we can write 

a{x,t) — a^(x,t) 8{x ~ mX)5{t — nT) 

m,n— — oo 

+a-{x,t) ^ {-lY5{x-raX)5{t-nT), (30) 



m,n— — oo 



where 

r+TT/X 

a^{x,t)= dfca±(fc)e*'="eT*"°('=)*. (31) 

J-tt/X 

Obviously, we can also write Ea. (|29|l in the form 

a[x, t) = a^[x, t) + (-1)" a^{x, t), (32) 

which will be used in the next section. 

Clearly the fields a^{x,t) and a'~{x,t) represent disturbances propagating 
(at least for small k) to the right and the left respectively, both with a dispersion 
relation coQ^k). The effect of the second dispersion relation Wi(fc) is to introduce 
a kind of "temporally antiferromagnetic" behavior, indicated by the (— 1)", 
associated with the leftward going wave. The solution seems the simplest when 
formulated in terms of these two fields a~^{x,t) and a~{x,t). Since we have 
dynamical equations, Eas. (|12(l . for both the Rm.n and the Lm,m we will have 
four fields in all as a ranges over R and L: i?+(x,t), R^{x,t), L^{x,t), and 
L-{x,t). 
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3.4 Initial conditions 

Eas. (|12|) requires two sets of initial conditions: the value of Rm,n and L^.n at 
both n = and n ~ 1. We can now see how these are built into the fields 
a'^^Xjt). We begin by noting that from Eas. l(T7|l and l|5n|l we can immediately 
identify 

a„j,o = a+ {mX, 0) + a" {mX, 0) , (33) 
ara.i = a+ {■mX,T) - a- {■mX,T) . (34) 

So to construct the a^(fc) appropriate for our initial conditions we must invert 
Eas. (|31|) to find the a^(fc) in terms of a^{x,t). This is easily done by focusing 
on points x = mX; writing Eas. (|31|l at such points and performing the sums 
indicated below, we find 

-J oo 

a± (fc) = — Xe=^''^°('=)* J2 a± (mX,i)e-*'=™^. (35) 

m— — OG 

As the left hand side is independent of time, so must be the right hand side. 
Writing down these equations for t = and t = T we can then use Eas. (|^ and 
l|34(l to solve for the a^(fc) in terms of the am,o and the am.i- The result is 



oo 



Aw ^ cos [Two 

m — — oo ^ 



47r v/l-psin^ kX 

where in the second expression we have used the dispersion relation (|26|) . 

3.5 Full solution 

The solution of the walk problem subject to arbitary initial conditions can then 
be written in the form of Ea. H30() or, identifying the discrete points of interest 
in space and time, in the form 

a™,„ = a+{mX,nT) + (-l)"a"(mX, nT), (37) 

where d^{x,t) are given by Eas. (|31|l with the Fourier components given by 
Ea. H36|l . Inserting those Fourier components we find that we can write the full 
solution for d^{x,t) in terms of a single Green function g{x;t), 

+ 00 +00 

a=^(x,t)= ^ g{±{x-mX);t-T)a,nfi± ^ g{±{x-mX);t)a„i^i, (38) 

m— — oo m— — oc 

where 

gi^^-t)^— dk =, (39) 

J-n/x ^/\- psm^kX 



8 



and we have used the fact that g{x;~t) = g{—x;t). Using the original walk 
equations (UJ and © we can then find the am.i in terms of the a„i_0: foi' a = R, L. 
Then our final solution is given by Eq.(|27I) with 

L=^(x,t) = ^g{±{x-mX);t-T)L,n,oTVp^g{±{x-mX);t)L,n+i.o 

7n m 

±^T~pY,9{H^~mX)-t)R„,^i^o , (40) 

and 

±^/T^Y.g{±ix - mX); t) Lm+i,Q ■ (41) 

m 

Earlier, Nayak and Vishwanath ^21 formally solved the discrete QW and 
analyzed the asymptotic behaviour of the solution. In the Appendix we show 
that our formal solution and theirs are the same. More recently, Carteret et al. 
[Tl] have reviewed previous results concerning the asymptotics of the QW and 
presented new asymptotics. Of course, all those asymptotic results fully apply 
to our solution; we do not pursue this further here. 

4 An approximate approach: the long wavelength 
limit 

The solution in the form we have constructed leads naturally to a passage to the 
continuum limit, which we identify by focusing on the continuous functions of 
space and time a{x,t), and constructing simplifed equations for their evolution 
in the long wavelength limit. 



4.1 Derivation of the equation of evolution in the long 
wavelength limit 

Our starting points are Eas. l|31|l and l|32(l . which we write again here for clarity: 

a{x,t) ^ a+{x,t) + {-!)" a- {x,t), (42) 
a^ix,t) = / dka^ {k)e'''''e^"^"^'''^\ (43) 

J-tt/X 

with a* (k) given by Eas. H36|l . We consider the long wavelength limit by ne- 
glecting high frequency (and small wavelength) effects. For that purpose we 
introduce low-frequency fields 

r+TT/L 

i± {x, t)= dkg (kX) a± (k) e^'^^e^F^^'^C^)*. (44) 

J-tt/L 
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which in the long wavelength limit take the place of the fields a^{x,t). Here 
Q{kX) is a function that vanishes at high spatial frequencies k; that is, it plays 
the role of a low-pass filter. The cut-off that it provides is assumed to be such 
that we can make a series expansion in luq (k). Taking the dispersion relation 
Eq. H26|l and expanding as a power series one gets 

LjT-^ij^T^ + ... = ^Xk-^^X^k"^ + ... (45) 

Using the lowest order, ujT — ^JpXk, we approximate (y>T)^ w (y^Xfc)^ and 
write to uj (fc) with 

Now, we take this aproximate dispersion relation in Ea. l44() and write 



^(fc)- VP^fc-^VP(l-p)^^'- (46) 



A^{x,t)K dka±{k)g{kX)e'''''e'^'^'^''^\ (47) 



where uq (fc) has been approximated by lo (fc) and the integration limits have 
been extended to infinity. Finally from Eas. H47(l it is easy to derive the equations 
of evolution of ix,t). The equations can be written in the form 



d I- p 



dS, 6 dC 



^±(C,r) (48) 



where we have introduced the new functions (^,r) = A^{X£,,Tt), and the 
normalized time and space variables r = t/T and ^ = x/X. Notice that the 
form of this equation is independent of the particular form chosen for the cut-off 
function Q (kX). 

This is the differential equation derived earlier (Eq.(15) of [H]) for the special 
case p = 1/2. It can be identified with the equation the envelope function of 
a light pulse would satisfy in a description of propagation through a material 
medium with no group velocity dispersion at the lowest order - such lowest 
order dispersion would appear through a term involving a second derivative 
with respect to ^ - but with higher order dispersion described by the third 
derivative term. The equation is well known in fiber optics, where it is used 
to describe pulse propagation near the zero-dispersion wavelength '16' . Notice 
that the effective "dispersion coefficient" in Eq. H48|l is proportional to (1 — p) 
and thus, interestingly, can be varied by modifying the unitary transformation 
applied in the QW after each displacement. 

4.2 Solution of the differential equation 

Let us proceed now to analyze the solution of Eq. H48|) . By Fourier transforming 
Ea. H48|l one easily gets 

/ + 00 
(q,0)e*«^e=F*''^^(l+^«')^ (49) 
-oo 
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with 



1 



(50) 



The initial conditions {^,0) can be calculated from Eas. ()36|l and (|47|l . 
First we rewrite Eas. (|35|l in the long-wavelength limit 



^ oo 

(fc) = — X {am,o ± ajn,i) e' 



2kX 



(51) 



where we have approximated exp [Tujq (k)] « 1. Substituting this into Ea. H47|l . 
and using the normalized spatial frequency q = Xk we get 



Jq{i-m) 



(52) 



m— — oo 



Now we need to choose a particular cut-off function. To obtain an explicit 
solution of Eq. we adopt a Gaussian, 



Q iq) - e- 



(53) 



where w is a free parameter. We further consider the usual case that the QW 
starts with the particle at position m = 0, and we get 



A±(^,0) = ao,oG(0)±a_i,iG(-l)±ai,iG(l), 



(54) 



G(m) 



N exp 



(55) 



with N a normalization factor that can be fixed by impossing 

for the total probability remains equal to unity at all times. In the following we 
omit this constant, as it is nothing but a scaling factor. 

In substituting the initial condition (|54|l into Ea. l5U() and this into Ea. H49|) . 
we are left with integrals of the form 

Z(^,r) = / dk eii-p[iAk + i-Bk^ -Ck' 



A 



1 



(56) 
(57) 
(58) 

C = u;^ (59) 
whose solution can be expressed in terms of the Airy function Ai (x) ^\ 

27r 



>3 = ^Vp{l-P)r, 



' ^ A,' 



\ 64/3 



(60) 
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Then the final result reads, up to a normalization factor, 



= (e,r) + (-!)" A- (C,r) (61) 
A±(C,r) = ao^o2(±^,T)±a_i^iZ(±(f + l),r)±ai,iZ(±(^-l),T), 



or, explicitely, 



with 



Ri^,r) = i?+(e,r) + (-l)"i?-(e,r) (62) 
Z(e,r) = i+(e,r) + (-l)"i-(e,r) (63) 



i?±(e,r) = i?o,o[2(±e,r)±V^Z(±(e-l),T)] (64) 

±v/r^Lo,o2 (± (e + l),T), 
L±(e,r) = io.o[2(±C,T)TVp2(±(e + l),T)] (65) 

±yWi?o,o2(±(C-l),r), 

where use was made of the discrete QW Eqs.(|7|l and 

Notice that while the form of Eq. H48|) is independent of the particular form 
of G (kX), this solution depends explictly on w because the application of the 
initial condition does depend on the details of that cut-off function. 

Finally, the probability distribution is given by 

P{^,r) = P''(i,T) + P^{^,T), (66) 
P^(e,T) - \A+itT) + i-irA-itr)\\ A^R,L. 

Let us now compare the above long wavelength approximation with the so- 
lution of the original discrete QW Eqs.Q and ((SJ. In the left column of Fig. 2 
we present the total probability distribution Pm (n) , as well as the partial 
probabilites P^ (n) and P^ (n), as obtained from Eqs.Q and © for n = 200, 
^0,0 = l/v^ and Lq^q = i/V^ (a value of the initial conditions for which the 
mean displacement is null) and p — 1/2 (which corresponds to the Hadamard 
walk). In the right column of Fig. 2 we present the solution of the long wave- 
length approximation P{^,t) (and P^{^,t) and P^ {^,t)) for t — 200 and 
the same value of p and of the initial conditions and for w = 0.4. Finally, in 
the central column of the figure we present the long wavelength approximation 
evaluated only at integer values of the position ^. These last plots are given 
in order to make a better comparison with the exact solution, as the latter is 
nonvanishing only at even (odd) positions at even (odd) times. Although we 
only show plots for a particular initial condition and value of p, similar corre- 
spondence between the exact solution and the long wavelength limit results for 
other sets of initial conditions and values of p as well. 

In Fig. 2 we chose the value w — 0.4 because for that particular value 
the agreement of the long wavelength approximation with the exact solution is 
close to the best that can be achieved. In Fig. 3 we plot the long wavelength 
approximation for three different values of w. Note that for the largest w — 0.55 
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the long wavelength approximation yields values that are too low, compared 
to the exact solution, for small while for the smallest w — 0.25 the long 
wavelength approximation extends to |^| beyond the range of the exact solution. 

We see from Fig. 2 that the main trends of the exact QW solution are well 
captured by the long wavelength approximation, although obvious quantitative 
differences appear. The main discrepancy between the exact solution and the 
long wavelength approximation is the presence of the oscillations that appear 
in the latter. By taking into account higher order terms in the expansion H45|l 
a more accurate approximation to the exact solution could be obtained, but 
the corresponding differential equation would not have a simple analytic solu- 
tion. In any case, our main thesis is not that Eqs. H62|l and H63|l are a good 
quantitative approximation to the exact solution. Rather, our central result is 
that comparisons such as those in Fig. 2 demonstrate that the propagation and 
dispersion effects described by Eq. H48(l capture the qualitative nature of the 
QW. 

5 Conclusions 

In this paper we have rigorously derived both the exact solution of the coined 
QW for arbitrary unitary transformations, and a long wavelength approximation 
that constitutes an approximate continuous limit of the QW. 

The long wavelength approximation identified here corresponds to the dy- 
namical equation of an envelope function characterizing a light pulse propa- 
gating in a medium with higher order dispersion, for example in an optical 
fiber close to the zero-dispersion wavelength |^. This, together with the fact 
that entirely classical implementations of the unidimensional QW are possible 
E|) lead us to the conclusion that entanglement does not play a central 
role in the unidimensional QW, and that the walk can be regarded as a pure 
wave phenomenon. This conclusion cannot, however, be generalized to higher 
dimensional QWs ^H]- There several coins are needed, and it would in such 
a scenario that the entanglement characteristic of an essential quantum nature 
might appear, as already discussed [10) . It would then be very interesting to 
generalize the work we have carried out here by deriving long wavelength de- 
scriptions of multidimensional QWs. This could help elucidate any role played 
by entanglement in such processes. 

Another point of interest is the effect of decoherence on QWs. Kendon and 
Tregenna [201 have shown that decoherence modifies the probability distribution 
of the QW in somewhat unsuspected ways. In particular, it can lead to quite 
smooth, and at the same time wide, probability distributions. It would be 
interesting to investigate the role of decoherence in continuous descriptions of 
the QW, such as that we have presented here. In this sense, the recent work 
by Romanelli et al. j21| . in which Markovian and intereference terms of the 
quantum evolution are separated, provides a new interesting tool for this type 
of analysis. 
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7 Appendix 

In this appendix we show the relation between the walk definition we use here 
and that used by Nayak and Vishwanath ll^l . and confirm that our solution 
reproduces theirs. We begin with our walk in the case of a Hadamard transfor- 
mation, 

Rm,n = (^m-l,ri-l + ^m+l,n-l) , (67) 

Lm,n — {Rm-l,n-l — im+l,ra-l) • (68) 

and define new variables 

-^m^n — — m+l,n; -^m.n — m— l,n: (69) 

so 

^m^n m+l,ni -^m.n — m— l,n. C^^) 

In terms of these new variables the walk Eas. (|67|) and (|68|l are specified by 

Lm,n = (j^m+l,n-l + ^m-|-l,n-l^ , (71) 

Rm,n = + Rm-l,n-l^ ■ (72) 

These are the walk equations of Nayak and Vishwanath |2| • Then, in terms of 
our functions a^{x,t), where a is i? or L, we have from Eas. (|37|l and (|69|l . 

RmJ^ = L+{{-m+l)X,nT) + {-iy'L-{{-m + l)X,nT), (73) 
W = i?+((-m-l)X,nT) + (-l)"i?-((-m-l)X,nr). (74) 

In writing out the expressions for and at the lattice sites indicated we 
do a change of variables, putting k = k' + tt/X. We use 

ujo{k' + Tr/X) = -ujo{k'), (75) 
exp[i{k' + tt/X)X] = -exp[ik'X], (76) 
exp[^i{k' + 7r/X)mX] = (-1)™ exp [-ifc'mX] , (77) 
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and can then adjust the hmits of integration of the function back to — tt/X to 
tt/X in the new variable k' , since the function is periodic over 2t: / X. Changing 
the name of the new variable back to k we can collect the terms and write 



Rn 



r/X 



dk e'^^ 



L+{k) - (_l)"+'"L-(fc + — ) 



-■k/X 
^^—ikmX^ — iujQ (k)nT 

dke^'^^ 

-tt/X 



X' 



i?+(fc)-(-l)"+™i?-(fc + -) 



^^—ikmX ^—iuiQ {k)nT 



(78) 



(79) 



Now the initial condition adopted in |13| is Lm,o = ^mo and Rm,o = 0, which, 
in our notation Eas. H70|) read Rm.o = ^m.-i and L^.o = 0. After one time step, 
see Eas. (|H7jl and we have Rm.i = Lm,i — Smo/V^, so we can construct 

L^{k) and i?=^(fc) from with p = 1/2. Wc find 



L+{k) 



-L-ik + -) 



X 



R^{k) = -i?-(fc + -) 



47r Vl + cos2 kX' 

X ^e^^o(k)T^tkx _,_ I 



4:n VI + cos2 kX 



and using H26|l we recover 

Rvn.Ai. — 



X 



1 + (-1)"-™J ^ 



tt/X 



JkX 



dk- 



tt/x Vl + cos^ kX 



^ — ikraX ^ — iujT){k)nT 



n+rn] 



X 
47r 



cos kX 



tt/X / 
-tt/X \ Vl+COS^fcX 



^—ikraX ^—iuq (k)nT 



in agreement with Nayak and Vishwanath 



(80) 
(81) 



(82) 



(83) 
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Figure Captions 



Fig.l. Dispersion relations as given by Ea. H26|l 

Fig. 2. Probability distribution of the QW. The exact solution is presented 
in the left column (P^ (n), P^^ (n) and P,f^ (n) for n = 200 from top to bottom), 
the long wavelength approximation (P (^, t), P^ (^, t), and {£,, t) for t = 200 
from top to bottom) is presented in the right column, and also in the central 
column but in this case evaluated only at integer even. In the left and central 
columns the points have been joined to guide the eye. The initial condition is 
-Ro.o = l/\/2 and io.o = i/V^, and the parameters are p = 1/2, (Hadamard 
walk) and w — 0.4 for the long wavelength approximation. The probabilities 
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in the long wavelength approximation have not been normalized, and the units 
are thus arbitrary. 

Fig. 3. P r) for r = 200 and the same initial conditions and parameters 
as in Fig. 3 except for the values of w that are the indicated in the figure. 
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